Kasteleyn cokernels 



Greg Kuperberg| 
UC Davitf] 

We consider Kasteleyn and Kasteleyn-Percus matrices, which arise in enumerating matchings of planar 
graphs, up to matrix operations on their rows and columns. If such a matrix is defined over a principal ideal do- 
main, this is equivalent to considering its Smith normal form or its cokernel. Many variations of the enumeration 
methods result in equivalent matrices. In particular, Gessel-Viennot matrices are equivalent to Kasteleyn-Percus 
matrices. 

We apply these ideas to plane partitions and related planar of tilings. We list a number of conjectures, sup- 
ported by experiments in Maple, about the forms of matrices associated to enumerations of plane partitions and 
other lozenge tilings of planar regions and their symmetry classes. We focus on the case where the enumera- 
tions are round or (/-round, and we conjecture that cokernels remain round or g-round for related "impossible 
enumerations" in which there are no tilings. Our conjectures provide a new view of the topic of enumerating 
symmetry classes of plane partitions and their generalizations. In particular we conjecture that a ^-specialization 
of a Jacobi-Trudi matrix has a Smith normal form. If so it could be an interesting structure associated to the 
corresponding irreducible representation of SL()i, C). Finally we find, with proof, the normal form of the matrix 
that appears in the enumeration of domino tilings of an Aztec diamond. 



1. INTRODUCTION 

The permanent-determinant and Hafnian-Pfaffian methods 
of Kasteleyn and Percus give determinant and Pfaffian expres- 
sions for the number of perfect matchings of a planar graph 
pl| . Although the methods originated in mathematical 
physics, they have enjoyed new attention in enumerative com- 
binatorics in the past ten years JToj , p"2| , p3[ |l6| , [34j, in particu- 
lar for enumerating lozenge and domino tilings of various re- 
gions in the plane. These successes suggest looking at further 
properties of the matrices that the methods produce beyond 
just their determinants or Pfaffians. 

In this article we investigate the cokernel, or equivalently 
the Smith normal form, of a Kasteleyn or Kasteleyn-Percus 
matrix M arising from a planar graph G. One theme of our 
general results in Sections 33 and [O] is that the cokernel is a 
canonical object that can be defined in several different ways. 
More generally for weighted enumerations we consider M up 
to the equivalence relation of general row and column opera- 
tions. If G has at least one matching, then t he se t of matchings 
is equinumerous with cokerM. (In Section 4.2 , we conjecture 



an interpretation of this fact in the spirit of a bijection.) The 
cokernel of M is also interesting even when the graph G has no 
matchings, a situation which we call an impossible enumera- 
tion. Propp proposed another invariant of M that generalizes 
to impossible enumerations and that was studied by Saldanha 
J22I , |25j], namely the spectrum of M*M. 

The idea of computing cokernels as a refinement of enumer- 
ation also arose in the context of Kirchoff 's determinant for- 
mula for the number of spanning trees of a connected graph. 
In this context the cokernels are called tree groups and they 
were proposed independently by Biggs, Lorenzini, and Mer- 
ris [§ |§ f§]. Indeed, Kenyon, Propp, and Wilson Q, gen- 



eralizing an idea due to Fisher [^J, found a bijection between 
spanning trees of a certain type of planar graph G and the 
perfect matchings of another planar graph G'. We conjecture 
that the tree group of G is isomorphic to the Kasteleyn-Percus 
cokernel of G'. 

In Section 5^ we study cokernels for the special case of 
enumeration of plane partitions in a box, as well as related 
lozenge tilings. We previously asked what is the cokernel of 
a Carlitz matrix, which is equivalent to the Kasteleyn-Percus 
matrix for plane partitions in a box with no symmetry imposed 
p2]]. We give two con jectures that together imply an answer. 
Finally in Section 5.3 we derive, with proof, the cokernel for 
the enumeration of domino tilings of an Aztec diamond. 
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2. PRELIMINARIES 

2.1. Graph conventions 

In general by a planar graph we mean a graph embedded in 
the sphere S 2 . We mark one point of S 2 outside of the graph as 
the infinite point; the face containing it is the infinite face. Our 
graphs may have both self-loops and multiple edges, although 
self-loops cannot participate in matchings. 
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2.2. Matrix algebra 

Let R be a commutative ring with unit. We consider ma- 
trices M over R, not necessarily square, up to three kinds of 
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equivalence: general row operations, 
M^AM 

with A invertible; general column operations, 
M^MA 

with A invertible; and stabilization, 
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M , 



Mi 



and its inverse. Any matrix M' which is equivalent to M under 
these operations is a stably equivalent form of M. 

A matrix A over R is alternating if it is antisymmetric and 
has null diagonal. (Antisymmetric implies alternating unless 
2 is a zero divisor in R.) We consider alternating matrices up 
to two kinds of equivalence: general symmetric operations, 

A i — ► BAB T 

with B invertible; and stabilization, 
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M , 



and its inverse. A matrix A' which is equivalent to A is also 
called a stably equivalent form of A. 

As a special case of these notions, elementary row opera- 
tion on a matrix M consists of either multiplying some row i 
by a unit in R, or adding some multiple of some row i to row 
j i. Elementary column operations are defined likewise. We 
define a pivot on a matrix M at the (i, j) position as subtracting 
M k j/M j j times row i from row k for all k ^ i, then subtracting 
M i k /M i j times column j from column k for all j ^ k. This op- 
eration is possible when M ; ■ divides every entry in the same 
row and column. In matrix notation, if M 1 1 = 1, the pivot at 
(1,1) looks like this: 



M = 



A deleted pivot consists of a pivot at (/, j) followed by deleting 
row i and column j from the matrix. The deleted pivot at ( 1 , 1 ) 
on our example M looks like this: 




M- 




M'-XY T . 



If A is an alternating matrix, we define an elementary sym- 
metric operation as an elementary row operation followed by 
the same operation in transpose on columns. We can similarly 
define a symmetric pivot and a symmetric deleted pivot. All of 
these operations are special cases of general symmetric matrix 
operations, and therefore preserve the alternating property. 

If R is a principal ideal domain (PID), then an n x k matrix 
M is equivalent to one called its Smith normal form and de- 
noted Sm(M). We define Sm(M) and prove its existence in 



Section ||. Note that if M' is a stabilization of M, then Sm(M') 
is a stabilization of Sm(M). 

If R is arbitrary, then we can interpret M as a homomor- 
phism from R k to R". In this interpretation M has a kernel 
kerM, an image imM, and a cokernel 

cokerM = R"/imM. 

If R is a PID, the cokernel carries the same information as the 
Smith normal form. Over a general ring R, only very special 
matrices admit a Smith normal form. Determining equiva- 
lence of those that do not is much more complicated than for 
those that do. In particular inequivalent matrices may have 
the same cokernel. However, over any ring R the cokernel is 
invariant under stable equivalence and it does determine the 
determinant detM up to a unit factor. A special motivation for 
considering cokernels appears when R = Z and M is square. 
In this case the absolute determinant (i.e., absolute value of 
the determinant) is the number of elements in the cokernel, 

| detM | = | cokerM |, 

when the cokernel is finite, while 

detM = 

if the cokernel is infinite. 

An alternating matrix A over a PID is also equivalent to 
its antisymmetric Smith normal form Sm„(A), which we also 
discuss in Section |[ Again cokerA determines Sm a (A). 

Remark. If M is a matrix over the polynomial ring F[x] over 
an algebraically closed field F, which is a PID, then the factor 
exhaustion method for computing detM Jl4| | actually com- 
putes the Smith normal form (or cokernel) of M. Thus the 
Smith normal form plays a hidden role in a computational 
method which is widely used in enumerative combinatorics. 

The basic version of the factor exhaustion method computes 
the rank of the reduction 

M®F[x]/(x-r) 

for all r £ F. These ranks determine detM up to a constant fac- 
tor if the Smith normal form of M is square free. It is tempt- 
ing to conclude that the factor exhaustion method "fails" if the 
Smith normal form is not square free. But sometimes one can 
compute the cokernel of 

M®F[x]/(x-r) k , 

for all r and k. This information determines cokerM, as well 
as detM up to a constant factor, regardless of its structure. 
Thus the factor exhaustion method always succeeds in princi- 
ple. 



3. COUNTING MATCHINGS 

Most of this section is a review of Reference [l6|. 
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3.1. Kasteleyn and Percus 



> ■ ■ < 



Let G be a connected finite graph. If we orient the edges of 
G, then we define the alternating adjacency matrix A of G by 
letting A { • be the number of edges from vertex i to vertex j 
minus the number of edges from vertex j to vertex i. If G is 
simple, then the Pfaffian Pf A has one non-zero term for every 
perfect matching of G, but in general the terms may not have 
the same sign. 

Theorem 1 (Kasteleyn). If G is a simple, planar graph, then 
it admits an orientation such that all terms in Pf A have the 
same sign, where A is the alternating adjacency matrix of G 

771; 



In general an orientation of G such that all terms in Pf A 
have the same sign is called a Pfaffian orientation of G. If an 
orientation of G is Pfaffian, then the absolute Pfaffian | Pf A 
is the number of perfect matchings of G. Kasteleyn's rule for 
a Pfaffian orientation is that an odd number of edges of each 
(finite) face of G should point clockwise. We call such an 
orientation Kasteleyn flat and the resulting matrix A a Kaste- 
leyn matrix for the graph G. Likewise an orientation may be 
Kasteleyn flat at a particular face if it satisfies Kasteleyn's rule 
at that face. Every planar graph has a Kasteleyn-flat orienta- 
tion, although it is only flat at the infinite face of G if G has an 
even number of vertices. Forming a Kasteleyn matrix to count 
matchings of a planar graph is also called the Hafnian-Pfaffian 
method [|15|. 

Percus [pi|] found a simplification of the Hafnian-Pfaffian 
method when G is bipartite. Suppose that G is a bipartite 
graph with the vertices colored black and white, and suppose 
that each edge has a sign + or — , interpreted as the weight 1 
or — 1 . Then we define the bipartite adjacency matrix M of G 
by letting M i • be the total weight of all edges from the black 
vertex i to the white vertex j. If G is simple, then the determi- 
nant detM has a non-zero term for each perfect matching of 
G, but in general with both signs. 

Theorem 2 (Percus). If G is a simple, planar, bipartite 
graph, then it admits a sign decoration such that all terms in 
detM have the same sign, where M is the bipartite adjacency 
matrix of G. 

In the rule given by Percus, the edges of each face of G 
should have an odd number of — signs if and only if the face 
has 4k sides. We call such a sign decoration of G Kasteleyn 
flat and the corresponding matrix M a Kasteleyn-Percus ma- 
trix. Every planar graph has a Kasteleyn-flat signing, although 
it is only flat at the infinite face if G has an even number of ver- 
tices. Forming a Kasteleyn-Percus matrix M is also called the 
permanent-determinant method. A Kasteleyn matrix A for a 
bipartite graph G can be viewed as two copies of a Kasteleyn- 
Percus matrix M: 



A = 






M 


—M 






If the graph G is not simple, then we may make it simple 
by tripling edges, as shown in Figure 111 The set of match- 
ings of the new graph G' is naturally bijective with the set of 



Figure 1: Tripling an edge in a graph. 



matchings of G. A more economical approach is to define 
a Kasteleyn matrix A or a Kasteleyn-Percus matrix M for G 
directly. In this case A t • is the number of edges from vertex 
i to vertex j minus the number from j to i, while M ; ■ is the 
number of positive edges minus the number of negative edges 
connecting ; and /'. 

A variant of the Hafnian-Pfaffian method applies to a pro- 
jectively planar graph G which is locally but not globally bi- 
partite. This means that G is embedded in the projective plane 
and that all faces have an even number of sides, but that G is 
not bipartite. An equivalent condition is that all contractible 
cycles in G have even length and all non-contractible cycles 
have odd length. 

Theorem 3. If a projectively planar graph G is locally but not 
globally bipartite, then it admits a Pfaffian orientation. 

The orientation constructed in the proof of Theorem [3] is 
one with the property that each face has an odd number of 
edges pointing in each direction. We call such an orientation 
Kasteleyn flat; it exists if G has an even number of vertices. 
(If G has an odd number of vertices, then every orientation 
is trivially Pfaffian.) We call the corresponding alternating 
adjacency matrix A the Kasteleyn matrix of G as usual. 

The constructions of this section, in particular Theorems [j] 
and ||, generalize to weighted enumerations of the match- 
ings of G, where each edge of G is assigned a weight and 
the weight of a matching is the product of the weights of its 
edges. We separately assign signs or orientations to G using 
the Kasteleyn rule (in the general case) or the Percus rule (in 
the bipartite case). The weighted alternating adjacency ma- 
trix A is called a Kasteleyn matrix of G. If G is bipartite, the 
weighted bipartite adjacency matrix M, with the weights mul- 
tiplied by the signs, is a Kasteleyn-Percus matrix of G. Then 
Pf A or detM is the total weight of all matchings in G. 



3.2. Polygamy and reflections 

We can use the Hafnian-Pfaffian method to count certain 
generalized matchings among the vertices of a planar graph G 
using an idea originally due to Fisher [[7]]. We arbitrarily di- 
vide the vertices of G into three types: Monogamous vertices, 
odd-polygamous vertices, and even-polygamous vertices. An 
odd-polygamous vertex is one that can be connected to any 
odd number of other vertices in a matching, while an even- 
polygamous vertex can be connected to any even number of 
other vertices (including none). 

If G is a graph with polygamous vertices, we can find a 
new graph G' such that the ordinary perfect matchings of G' 
are bijective with the generalized matchings of G. The graph 
G' defined from G using a series of local moves that are shown 
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Figure 2: Resolving polygamy into monogamy. 



in Figure [2| (In this figure and later, an open circle is an even- 
polygamous vertex and a dotted circle is an odd-polygamous 
vertex.) We also describe the moves in words. First, if a 
polygamous vertex of G has valence greater than 3, we can 
split it into two polygamous vertices of lower valence with 
the same total parity. This leaves polygamous vertices of both 
parities of valence 1, 2, and 3. If a polygamous vertex v is 
even and has valence 1, we can delete it. If it is odd and has 
valence 1 or 2, it is the same as an ordinary monogamous ver- 
tex. If it is even and has valence 2, we can replace it with 
two monogamous vertices. If it is even and has valence 3, 
we can split it into an odd-polygamous divalent vertex and an 
odd-polygamous trivalent vertex. Finally, if it odd and has va- 
lence 3, we can replace it with a triangle. Each of these moves 
comes with an obvious bijection between the matchings be- 
fore and after. Thus these moves establish the following: 

Proposition 4 (Fisher). Given a graph G with odd- and even- 
polygamous vertices, the polygamous vertices can be replaced 
by monogamous subgraphs so that the matchings of the new 
graph G' are bijective with those of G. If G is planar, then G' 
can be planar. 

We call the resulting graph G' a monogamous resolution of 
G. If G is planar, then G' admits Kasteleyn matrices, and we 
call any such matrix a Kasteleyn matrix of G as well. 

■X- — > ■ < 




Figure 3: Moves on monogamous resolutions of a polygamous 
graph. 



The monogamous resolution of a polygamous graph is far 
from unique. But we can consider moves that connect dif- 
ferent monogamous resolutions of a polygamous graph. The 
moves are as shown in Figure [|: Doubly splitting a vertex, 
rotating a pair of triangles, and switching a triangle with an 
edge. Each of these moves comes with a bijection between 
the matchings of the two graphs that it connects. 

Proposition 5. Any two monogamous resolutions of a graph 
G are connected by the moves of vertex splitting and its in- 
verse, switching triangles, and switching a triangle with an 
edge. The moves also connect any two planar resolutions of a 
planar graph G through intermediate planar resolutions. 

The proof of Proposition || is routine. 




Figure 4: Removing a self-connected triangle. 

Another interesting move is removing a self-connected tri- 
angle, as shown in Figure [|. This move induces a 2-to-l map 
on the set of matchings before and after. 

Polygamous matchings have two common applications. If a 
graph G is entirely polygamous, then we can denote the pres- 
ence or absence of each edge by an element of Z/2. Each 
vertex then imposes a linear constraint on the variables, so the 
number of matchings is therefore either or 2" for some n. 
The corresponding weighted enumerations are related to the 
Ising model |Q, |l6| ^3|] . Another way to see that the number 
of matchings is a power of two is to use the moves in Figures || 
and^ to reduce a monogamous resolution of G to a tree, which 
has at most one matching. 




Figure 5: Using polygamy to count reflection-invariant matchings. 

Another application is counting matchings invariant under 
reflections [ |l5[ Suppose that a planar graph G has a re- 
flection symmetry a, and suppose that the line of reflection 
bisects some of the edges of G. Then the a -invariant match- 
ings of G are bijective with a modified quotient graph G//<7 
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in which the bisected edges are tied to a polygamous vertex, 
as in Figure ||. The parity of the polygamous vertex should be 
set so that the total parity (odd-polygamous plus monogamous 
vertices) is even. The same construction works if we divide G 
by any group acting on the sphere that includes reflections, 
since all of the reflective boundary can be reached by a single 
polygamous vertex. 



3.3. Gessel-Viennot 

The Gessel-Viennot method [|| ^] yields another determi- 
nant expression for a certain sum over the sets of disjoint paths 
in an acyclic, directed graph G. (Theorem ^ below, which is 
the basic result of the method, was independently found by 
Lindstrom [ 17 1- Gessel and Viennot were the first to use it 
for unweighted enumeration.) The graph G need not be pla- 
nar. We label some of the vertices of G as left endpoints and 
some as right endpoints, and we separately order the left end- 
points and the right endpoints. Let P be the set of collections 
of vertex-disjoint paths in G connecting the left endpoints to 
the right endpoints. If P is non-empty then there are the same 
number of left and right endpoints on left and right; if there 
are n of each we call the elements of P disjoint n-paths. The 
Gessel-Viennot matrix V is defined by setting V ( ■ to the num- 
ber of paths in G from left endpoint i to right endpoint j. 



Theorem 6 (Lindstrom, Gessel-Viennot). Let G be a di- 
rected, acyclic, weighted graph with n ordered left endpoints 
and n ordered right endpoints. Let P be the set of disjoint n- 
paths in G connecting left to right. IfV is the Gessel-Viennot 
matrix of G, then 



detV = JXO(-l)* 



(1) 



Here (— 1) is the sign of the bijection from the left to the right 
endpoints induced by the paths in the collection I, and w{€) is 
the product of the weights of the edges ofG that appear in £. 

Proof. We outline a non-traditional proof that will be useful 
later. We first suppose that the left endpoints are the sources 
in G (the vertices with in-degree 0) and the right endpoints 
are the sinks (the vertices with out-degree 0). We argue by 
induction on the number of transit vertices, meaning vertices 
that are neither sources nor sinks. 





Figure 6: Splitting a transit vertex. 



If G has no transit vertices, every path in G has length 
one. Consequently the n-paths in G are perfect matchings, 
and equation (Jl]) is equivalent to the definition of the determi- 
nant. Suppose then that p is a transit vertex in G. We form a 
new graph G' by splitting p into two vertices q and r, with q a 



sink and r a source, as shown in Figure ^. We number q and 
r as the n + 1st (last) source and sink in G' . We give the new 
edge between q and r a weight of — 1 . There is a natural bijec- 
tion between disjoint n-paths £ in G and disjoint n + 1 -paths 
£' in G'\ Every path in I which avoids p is included in £'. If 
some path in £ meets p, we break it into two paths ending at q 
and starting again at r. If £ is disjoint from p, we include the 
edge from r to q in £'. In order to argue that the right side of 
equation (jl]) are the same for G and G', we check that 

{-\) ( w{l) = (-\) ( 'w{£'). 

If £ avoids p, the two sides are immediately the same. If £ 
meets p, then (— 1) and (— 1) have opposite sign and so do 
w{£) and w{£'). The left side of equation ([l]) is also the same: 
If V and V ' are the Gessel-Viennot matrices of G and G', V is 
obtained from V' by a deleted pivot at (n + 1 , n + 1 ). 

Now suppose that the left and right endpoints do not co- 
incide with the sources and sinks. If G has a left endpoint q 
which is not a source, then there is an edge e from a vertex 
p to the vertex q. Let G' be G with e removed and let V' be 
its Gessel-Viennot matrix. Since the edge e is not in any n- 
path in G, the graph G' has the same n-paths with the same 
weights. If p is the ith left endpoint, we can obtain V' from V 
by subtracting w(i,j) times row i from j, where w(i,j) is the 
total weight in G of all paths from i to j. These row operations 
do not change the determinant. The same argument applies if 
G has a right endpoint which is not a sink. 

Finally if G has a right endpoint source or a left endpoint 
sink, then it has no n-paths and some row or column of V is 0. 
If G has a source or a sink which is not an endpoint, we can 
delete it without changing the Gessel-Viennot matrix V or the 
set of n-paths. □ 

We call the graph G' constructed in our proof of Theorem || 
the transit-free resolution of G. The transit-free resolution 
is a connection between the Gessel-Viennot method and the 
permanent-determinant method: 

Corollary 7. Let G be a connected, planar, directed, acyclic 
graph with n left and right endpoints on the outside face. Sup- 
pose that the left endpoints are segregated from the right end- 
points on this face. If the left endpoints are the sources and the 
right endpoints are the sinks, then the Gessel-Viennot matrix 

V of G is obtained from a Kasteleyn-Percus matrix M of the 
transit-free resolution G' of G by deleted pivots. If not every 
left endpoint is a source or not every right endpoint is a sink, 

V is obtained by deleted pivots and other matrix operations. 

Note that by construction the matchings of the transit-free 
resolution G 1 of G are bijective with the n-paths in G. The 
planarity of G together with the position of its endpoints imply 
that all n-paths induce the same bijection and therefore have 
the same sign. 

Proof. The Gessel-Viennot matrix V' of the graph G' is ob- 
tained from V by the stated operations. Thus it suffices to 
show that G' is planar and that V' is also a Kasteleyn-Percus 
matrix of G'. 
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Figure 7: Left-to-right orientation implied by segregation of sources 
and sinks (circled). 



We first establish that the orientation of G is qualitatively 
like that of the example in Figure ^: the orientations all point 
from left to right. More precisely, the edges incident to each 
transit vertex are segregated, in the sense that all incoming 
edges are adjacent and all outgoing edges are adjacent. The 
edges of each internal face are also segregated, in the sense 
that the clockwise edges are adjacent and the counterclock- 
wise edges are adjacent. To prove that G has this structure, we 
reallocate the Euler characteristic of the sphere, 2, expressed 
as a sum over elements of G. In this sum, each vertex and face 
has Euler characteristic 1 and each edge has Euler character- 
istic — 1 . If a pair of edges shares both a vertex v and a face 
/, we deduct i from the Euler characteristic of / if the edges 
both point to or both point from v, and otherwise we deduct 
j from v. Since each edge participates in 4 such pairs, these 
deductions absorb the total Euler characteristic of all edges. 

The reallocated characteristic of a vertex is 1 if it is a source 
or sink, if it is a segregated transit vertex, and negative other- 
wise. The reallocated characteristic of a face is at most 2 — 2n 
if it is the outside face (since orientations must switch between 
clockwise and counterclockwise both at the sources and sinks 
and between them), if it is a segregated internal face, and 
negative otherwise. (No face has positive reallocated charac- 
teristic since G is acyclic.) Thus the only way that the total 
can be 2 is if all internal faces and all transit vertices are seg- 
regated. 

That the transit vertices of G are segregated implies that G' 
is planar. That each internal face / is segregated implies that 
if / has k sides, the corresponding face /' of G' has 2k ~ 2 
sides. Moreover the k — 2 new edges of /' have weight —1 
in the proof of Theorem ^, which agrees with the Kasteleyn- 
Percus rule. Thus V' is a Kasteleyn-Percus matrix of G', as 
desired. □ 

Finally we have not discussed a Pfaffian version of the 
Gessel-Viennot method defined by Stembridge [29). We be- 
lieve that this method can be generalized further, and that it 
admits an analogue of Corollary [7| 



4. MATCHINGS AND SMITH NORMAL FORM 

4.1. Equivalences of Kasteleyn and Kasteleyn-Percus matrices 

If M is a Kasteleyn-Percus matrix of a bipartite, planar 
graph G, then we can consider its cokernel, which by Sec- 



tion 2.2 is equinumerous with the number of matchings of M 
if it has any matchings. Furthermore, if cokerM is infinite or 
if M isn't square, we can think of cokerM as a way to "count" 
matchings in a graph that has none. We call such a compu- 
tation an impossible enumeration. Both observations are rea- 
sons to study cokerM as part of enumerative combinatorics. 

If G is weighted by elements of some ring R, then we can 
consider M up to stable equivalence, whether or not it has a 
Smith normal form. 

Suppose that M and M' are two Kasteleyn-Percus matrices 
for the same planar graph G. Then the signs on G given by M' 
and M differ by a 1-cocycle c with coefficients in the group 
{+, — } Jig], Since the sphere has no first homology, c = 8d, 
where d is a 0-cochain. More explicitly d is a function from 
the vertices of G to {+, — }. We can use d to form two diag- 
onal matrices A and B with diagonal entries ±1 and such that 
M' = AMB. Evidently A and B are invertible over Z, so M' 
and M have the same cokernel. In conclusion: 

Proposition 8. If G be a weighted bipartite planar graph, 
then all of its Kasteleyn-Percus matrices M are stably equiva- 
lent forms. In particular cokerM is an invariant of G. 





Z/30Z/3 





zez/3 



Figure 8: Embedding-dependent Kasteleyn-Percus cokernels. 



Example 9. The Smith normal form or cokernel ofM can de- 
pend on the embedding of G in the plane. The top two graphs 
in Figure ^have cokernels Z/9 and Z/3 ©Z/3 in the two em- 
beddings shown. If G has an odd number of vertices, then 
it cannot be Kasteleyn flat on its outside face. In this case 
changing which face is on the outside can change the cokernel 
as well. The bottom two graphs in Figure [| are an example. 

Our analysis generalizes to the non-bipartite case. If G is a 
planar graph with a Kasteleyn matrix A, then we can consider 
A up to equivalence. 

Again all Kasteleyn matrices we choose for the planar 
graph G are equivalent, because any two Kasteleyn-flat ori- 
entations of G differ by the coboundary of a 0-cochain on G 
with values in {+, — }. The matrices are consequently equiv- 
alent under the transformation 

A i ► B T AB 
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for some diagonal matrix B whose non-zero entries are ±1. 
We can also pass from the usual clockwise-odd Kasteleyn rule 
to the counterclockwise-odd rule by negating A. We have no 
reason to believe that A and —A are equivalent over a general 
ground ring R, but they do have the same cokernel and are 
therefore equivalent if R is a PID. 

If G is projectively planar and locally but not globally bi- 
partite, the argument is slightly different while the conclusion 
is the same. In this case the cohomology group 



; ,Z/2) ^Z/2 



is non-trivial. Suppose that we have two Kasteleyn-flat orien- 
tations of G whose matrices are A and A'. Their discrepancy 
is a 1-cocycle c which could represent either the trivial or the 
non-trivial class in H 1 (SUP 2 ,Z,/2). If c is trivial, then 

A' = B T AB 

for some diagonal B. I.e., A and A' are equivalent. If c is 
non-trivial, then 

A' = —B T AB, 
i.e., A' is equivalent to —A. 



-X- 



Figure 9: Preserving Kasteleyn flatness. 



Kasteleyn and Kasteleyn-Percus matrices remain equiva- 
lent under more operations than just the choices of signs or 
orientations. In part icular they remain equivalent under the 
moves in Section 32. In each move we make a graph G' from 
the graph G, and we need to choose related Kasteleyn-flat ori- 
entations of both graphs. For example consider a double ver- 
tex splitting. If G is Kasteleyn-flat, and if we orient the two 
new edges in the splitting in opposite directions as in Figure |J 
then G' is also Kasteleyn flat. If the three vertices are num- 
bered 1, 2, and 3, then the matrix A' of G' has a submatrix of 
the form 




If we perform a deleted pivot at (1,2) and (2, 1 ), we reduce A' 
to the matrix A of G. 



4.2. Is it bijective? 

Whenever two sets are known to have the same size, a tra- 
ditional question in combinatorics is whether or not there is a 
bijection between them. In this section we conjecture a rela- 
tionship between cokernels and matching sets which is similar 
to a bijection. 



If M is a Kasteleyn-Percus matrix for an unweighted bipar- 
tite, planar graph G with at least one perfect matching, then 
two such sets to consider are cokerM and P, the set of per- 
fect matchings of G. To this pair we must add a third set, 
cokerM 7 , since the choice between M and M T is arbitrary. 
As explained in Section g cokerM and cokerM 7 are isomor- 
phic, but there is no canonical isomorphism. This is evidence 
against a natural bijection between cokerM and cokerM 7 , and 
therefore a natural bijection between either of them and P. On 
the other hand, the special planar structure of M might yield 
such bijections. 

It may be better to consider quantum bijections or linearized 
bijections instead of traditional ones. If A and B are two finite 
sets, a quantum bijection is a unitary isomorphism 

C[A] = C[B] 

between the formal linear spans of A and B. A quantum bijec- 
tion can be implemented by a quantum computer algorithm 
just as a traditional bijection can be implemented on a stan- 
dard computer p8[ | . A linear bijection is a linear isomorphism 

F[A] = ¥[B], 

not necessarily unitary, for some field F. A linear bijection 
does not have the empirical computational interpretation that 
a traditional bijection or a quantum bijection does, but as a 
means of proving that A and B are equinumerous, it can be 
considered constructive. 

If M is a non-singular n x n matrix over Z, then there is 
a natural quantum bijection between cokerM and cokerM 7 , 
namely the discrete Fourier transform. (It is also a special 
case of Pontryagin duality [^3|]). We express it by defining a 
unitary matrix U whose rows are indexed by x G cokerM and 
whose columns are indexed by y G cokerM 7 . Given such x 
and y, we let X and Y be lifts in Z" . We then define 



U X ,y = 



y^fdetMl 



It is easy to check that Y T M~ l X changes by an integer if we 
change the lift X of x, because two such lifts differ by an ele- 
ment in imM. 

Given M, G, and P as above, it might be possible to fac- 
tor the unitary map U into maps to and from C[P]. However 
we may need to further relax the notion of a bijection. Some- 
times when G is a finite group equinumerous with a finite set 
S, there is no natural bijection between them, but instead there 
is a natural, freely transitive group action of G on S. Hav- 
ing introduced quantum bijections, we can try to make C[P] 
a free unitary module over the group algebras C [cokerM] and 
C [cokerM 7 ]. We can even ask that the two group actions be 
compatible with U by requiring the commutation relations 

a y a x = exp(2niY T M~ l X)a x a y , 

where x G cokerM and y G cokerM 7 and a x and a y are their 
hypothetical actions on C[P], A standard theorem in represen- 
tation theory says that for any M the algebra 

D = C [cokerM] <g> C [cokerM 7 ] 
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twisted by this commutation relation is isomorphic to a matrix 
algebra. This means that D has only one irreducible represen- 
tation, and we conjecture that C[P] has the structure of this 
representation. 

If A is a non-singular alternating matrix, then we can de- 
fine a similar algebra D as a deformation of the group algebra 
C [coker A] . We let D be the formal complex span of elements 
a x with x £ coker A, and we arbitrarily order the elements in 
coker A. For x < y 6 coker A, we define 

OCxOCy — CX x +y, 

and for arbitrary x and y, we impose the relation 

<Xy<x x = exp(2niY T A~ l X)a x a^,. 

The algebra D is again isomorphic to a matrix algebra, since 
Theorem [l^ allows us to put A into 






M ' 


—M 






If we do so then algebra D then has the form previously de- 
scribed. 

Conjecture 10. If A is a Kasteleyn matrix of a planar graph 
G with non-empty matching set P, then there is a natural ac- 
tion of the algebra D on C[P], possibly depending on the way 
that G is embedded in the plane. 

5. LOZENGE AND DOMINO TILINGS 
5.1. Plane partitions and lozenge tilings 

Plane partitions are an interesting source of enumerative 
planar matching problems. We consider the cokernels and 
Smith normal forms that arise in these problems, relying on 
the material in References |l5| and [To] . 




Figure 10: A plane partition and a tiling. 

A plane partition in an a x b x c box is equivalent to a 
lozenge tiling of an (a,Z?,c)-semiregular hexagon (Figure [ic|), 
which in turn is equivalent to a perfect matching of the dual 
hexagonal graph Z(a,b,c) (Figure |TT| ). 

If a group G acts on the box and the plane partitions inside 
it, it also acts on the graph Z(a,b,c). The G-invariant match- 
ings correspond to matchings of a modified quotient graph 
Z G (a,b,c). To understand these graphs we recall the three 
generators of G: 



Figure 1 1 : The graph Z(2, 2, 3) . 



• p, cyclic symmetry for plane partitions or rotation by 120 

degrees for lozenge tilings, defined when a = b = c. 

• T, symmetry for plane partitions or diagonal reflection for 

lozenge tilings, defined when b = c. 

• K, complementation for plane partitions or rotation by 180 

degrees for lozenge tilings. 

We describe Z G (a,b,c) case by case: 

• If G = (p) or (p,k), or if a, b, and c are all even and 

G is (k), then Z G (a,b,c) is the usual quotient graph 
Z(a,b,c)/G. 

• If G — (kt) or G = (KT,p), then we delete the edges and 

vertices of Z(a,b,c) along the lines of reflection and let 
Z G (a,b,c) be a connected component of the remainder. 

• If G = (k) and only one or two of a, b, and c is even, then 

Z(a,b,c) has a central edge e invariant under K. If one 
dimension is even, we define Z K (a,b,c) by deleting e 
and its vertices and then quotienting by k. If two di- 
mensions are even, we define Z K (a,b,c) by deleting e 
but not its vertices, and then quotienting by K. 

• If G = (t), (t, k), or (t,P,k), we start with Z H (a,b,c), 

where H is, respectively, (1), (tk), or (p,Tff). We 
cut Z H (a,b,c) by the line of reflection, and tie the cut 
edges of one region to a polygamous vertex to define 
Z G (a,b,c). 

• If G = (T,p), we cut Z(a,a,a) by three lines of reflec- 

tion and tie the cut edges of one region (lying along 
two of the three lines) to a polygamous vertex to form 
Z G (a,a,a). 

For every symmetry group G, the number of G-invariant 
plane partitions N G (a,b,c) is round, meaning a product of 
small factors. (This is not a completely rigorous notion. In 
this case N G (a,b,c) grows exponentially in ab + ac + be while 
its prime factors grow linearly in a + b + c.) In addition if G 
is a subgroup of (p, t), the G-invariant partitions have round 
^-enumerations, where the g-weight of a plane partition is 
sometimes the number of cubes and sometimes the number 
of G-orbits of cubes. (To say that a polynomial P(q) is round 




Figure 12: The graph Z(2, 2, 3) weighted for ^-enumeration. 



means not only that it is a product of small factors, but also 
that the factors are cyclotomic. Equivalently P(q) is a ratio of 
products of differences of monomials.) All of these enumera- 
tions are proven [Q, JO], [H]] except for the conjectured or- 
bit ^-enumeration of totally symmetric plane partitions. Both 
orbit and cube ^-enumerations can be realized by weighted 
enumerations of matchings in Z G (a,b,c) jl5| , [lq ]. Figure |l2| 
shows an example where G is trivial (so that there is no dis- 
tinction between cube and orbit ^-enumeration in this case). 
Also the weight of a matching in the example only agrees 
with the ^-weight of the corresponding plane partition up to 
a (matching-independent) factor of q. We will therefore con- 
sider ^-enumerations over the ring 7L\q,q~ x \ and absorb pow- 
ers of q in normalization. 

We let Z G (a,b,c;q) be Z G (a,b,c) with weights chosen for 
cube ^-enumeration, and we let Z G (a,b,c;q) be Z G (a,b,c) 
with weights chosen for orbit ^-enumeration. 




A 



Figure 13: A lozenge tiling of a region with a = 4, A = (2,2), and 
M = (l)- 



The problem of counting matchings in Z(a,b,c) has two 
other interesting generalizations. The lozenge tilings of a par- 
allelogram strip with notches on both sides (see Figure |l3|) 
are naturally bijective with the semi-standard skew tableaux 
of shape Xj\i and parts bounded by some a §7.10], 
where X and ji are two partitions with X containing ji. If 
X has b nonzero parts, the parallelogram then has dimen- 
sions a by Aj + b. On the bottom row it has notches at po- 
sitions Xj + b+l — i, counting from the left. On the top row 
it has notches at positions fi, + b + 1 — ;, counting from the 
left and extending ji by so that it also has b parts. We leave 
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Figure 14: Weights of Z(X/fi;x). 



the proof of the bijection between tilings and tableaux to the 
reader since it is similar to existing arguments in the literature. 
We let Z(X/fi,a) be the graph dual to the tiling of this re- 
gion by triangles, so that the lozenge tilings correspond to the 
matchings of Z(X/ji,a). In an important weighted enumera- 
tion we assign the weight x. to the northeast-pointing edges in 
the rth row, as in Figure [l4| Call the weighted graph Z(X / ' ji;x) 
The total weight of its matchings is the skew Schur function 
s^ifS?) ( m finitely many variables). If we let 

q a = (l,q,... ,q"~ 1 ), 

then the specialization ^^(^a) is the standard ^-enumeration 
of skew tableaux. In particular, if we omit jj, (by setting it 
to the empty partition), then (x) is the character of an ir- 
reducible representation V(X) of GL(a,C), while s«(q a ) is 
round by the g-Weyl dimension formula. 

A second generalization is to count tilings of a semireg- 
ular hexagon with side lengths a,b + d,c,a + d,b,c + d and 
with a triangle of size d removed We do not know of 
a ^-enumeration of these, although the symmetry fCT (corre- 
sponding to TCPPs) appears when b = c and the symmetry 
p (corresponding to CSPPs) appears when a = b = c. We let 
Z(a,b,c,d,d) be the dual graph and we analogously define 
Z G (a,b,c,d,d) if G is a subgroup of (jCT,p). 



a a + d 




a + d 



Figure 15: Hexagonal regions without central triangles. 

As we mentioned in Section |], we can consider the Kaste- 
leyn or Kasteleyn-Percus cokernel for impossible enumera- 
tions where there are no matchings. This allows us to vary 
the above graphs in several ways, which we describe case by 
case: 
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• We can consider Z K (a,b,c) when all three of a, b, and c are 

odd, as well as Z^ K p Ja,a,a) when a is odd. 

• If one of a, b, and c is odd, we define the graph Z' K (a,b, c) by 

removing the central edge e but not its vertices before 
quotienting by K. If two are odd we define Z' K (a,b,c) 
by removing e and its vertices before quotienting by K. 

• If Z G (a,b,c) has a polygamous vertex, we can give it 

the wrong parity to make Z' G (a,b,c). If G = (t) or 
G = (p,T), then we also define the g-weighted forms 
Z' G (a,b,c;q) and Z' G (a,b,c;q). 

• We define the graph Z(a,b, c, d, e) by removing a triangle of 

size e from a semiregular hexagon of side lengths a, b + 
d,c,a + d,b,c + d, with d ^ e, as shown in Figure |l5[ 
The parameters d and e might even have opposite sign, 
which we interpret by turning the triangle upside-down, 
as also indicated in Figure |l5[ 

For each of the graphs defined in this section, we denote 
a corresponding Kasteleyn-Percus matrix by replacing Z by 
M if it is bipartite and monogamous, and the corresponding 
Kasteleyn matrix by replacing Z by A otherwise. 

Conjecture 11. Each of the matrices M(a,b,c;q), 
Mp(a,a,a;q), A x (a,b,b;q), A x (a,b,b;q), A' z (a,a,b;q), 
A' T (a,a,b;q), A^ p ^(a,a,a;q), A' <p ^(a,a,a;q), and M(p;q a ) 

admits a Smith normal form over Jj[q,q~ l ], and the entries 
are q-round. The Smith normal form over 1 of each of the 
matrices M(a,b,c,d,e), M KZ (2a,b,b), M^ p K ^ (2a, 2a, 2a), 

A K (a,b,c), A' K (a,b,c), A {Kt) (2a,b,b), A' {Kz) (2a,b,b), 

has round entries. 

We could have stated conjecture [Tl] in greater generality 
by combining more of the variations above. For example the 
graph Z(a,a,a,d,e) has symmetries and we could consider 
the matrix of a suitably modified quotient graph. Indeed we do 
not know how to state Conjecture |ll] in full generality, since 
there are yet other variations of counting lozenge tilings in a 
hexagon with round enumerations [J4J]. Presumably many or 
all of these variations also have impossible counterparts, and 
some may have ^-enumerations as well. Also the spirit of the 
conjecture is to find the Smith normal forms or the cokernels 
explicitly. 

Conjecture 12. The Smith normal forms over 1\q,q~ l \ of 
each of the matrices M(a,b,c;q), Mp(a,a,a;q), A z (a,b,b;q), 
A T (a,b,b;q), andA^ p ^(a,a,a;q) are square free. 

Conjecture [l2| may also not be fully general, although we 
note that the Smith normal form of M(X;3a) is not always 
square free. Note that Conjectures |l2| andO would together 
solve problem 5 in Propp's problem list [|22|], which asks for 
the Kasteleyn cokernel of the graph Z(a,b,c). The reason is 
that the Smith normal form of any non-singular matrix M is 
uniquely determined by detM, provided that the normal form 



exists and is square free. For example the ^-enumeration of 
plane partitions tells us that 

detM(2,2,2;<7) = (2) 2 q (5) q , 

where 



So Conjectures |TT] and[j~2| assert that Sm(M(2,2,2;g)) exists 
and its nontrivial entries are (2) q and (2) q (5) q . If we special- 
ize at q = 1 , we obtain the correct prediction that 

cokerM(2,2,2) ^Z/2®Z/10. 

Stembridge [ ^o| ] noticed a family of relations, called the 
q = — 1 phenomenon, between — 1 -enumeration of some sym- 
metry classes of plane partitions and ordinary enumeration of 
other symmetry classes. We conjecture an extension of the 
phenomenon that is easier to state in terms of cokernels in 
some cases and Smith normal forms in others. 

Conjecture 13. If G is (1) or (p) and ifG' = (G, k), then 

Sm(A G ,(a,b,c)) = Sm(M G (a,b,c)_ l ) 2 . 

IfG is (t) or (p, t), and if G 1 is, respectively, (kx) or (p, Kl), 
then 

cokerA G (a,fo,c)_j = cokerM G ,(a,b,c)® 2 . 
IfG is (t) or (p, t) and ifG' — (G, k), then 

cokerA G ,(a,b,c) = cokerA' G (a,b ,c) _ Y . 

5.2. Jacobi-Trudi matrices 

One conclusion of our construction for skew tableaux is a 
novel determinant formula for skew Schur functions: 

s X/n@) =detM(X/n;x). 

We can compare this to two other determinant formulas for 
skew Schur functions [ |27| , §7.16]. Let 

J(X/li;Z) iJ = h Xi _ flj _ i+ .(x) ) 

where h n (x) is the «th complete symmetric function of x, and 
let 

D(X/ii;?) iJ = e xh ^_. +j .(x), 

where e„ (x) is the «th elementary symmetric function of x and 
A' is the partition conjugate to X. Then the Jacobi-Trudi iden- 
tity states that 

S XI^) =det7(A/^i;x) 
while the dual Jacobi-Trudi identity states that 
s x/^x) =detD(X/n;x). 
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Figure 16: First Gessel-Viennot model of lozenge tilings. (Not all 
orientations are shown.) 



one can check that the transit-free resolution of Y(X /jU;3t) is 
Z(X / jx;x). Thus Corollary says that J(X / jx;x) is equivalent 
toD(X/fi;x). □ 

Question 15. Are J{X / H',x) and D(X / 1 li;x) stably equivalent 
over the ring of symmetric functions? 

Theorem |l4| implies that there are several ways that equiv- 
alent forms of the matrix M(X/fi;x) arise in several com- 
mon guises. If we set x = q a and \i = 0, then Conjec- 
ture [ll] asserts that M(X;q a ) admits a Smith normal form a 
Smith normal form over the ring l\q,q~ ]. This suggests that 
cokerM(A;g„) is an important extra structure that one can 
associate to the representation V(X) mentioned previously. 
Moreover the relationship between M(X) and V (X) is an in- 
teresting special case of Conjecture |To[ 



Theorem 14. The Jacobi-Trudi matrices J(X/^i;x) and 
D(X /ju;3c) are stably equivalent over the ring Z[x] to the ma- 
trix M(X I fi;x). 

Proof. The theorem is a special case of Corollary ^. Let 
X(X/jJ.;x) be the graph exemplified in Figure |l6| a square 
grid of height a — 1 and width X 1 + b. On the top there is a left 
endpoint at each position [L + b + 1 — i and on the bottom a 
right endpoint at each position A- + b + 1 — i. Each horizontal 
edge points to the right and on the z'th row has weight x ( . Each 
vertical edge points down. By a standard argument the Gessel- 
Viennot matrix of X(X / jx;x) is J(X /n;x). At the same time, 
one can check that the transit-free resolution of X(X / 'jJ,;x) is 
Z(X /fi;x). Thus Corollary [7] says that J(X /jx;x) is equivalent 
toM(X/n;x). 



5.3. Domino tilings 




Figure 17: Second Gessel-Viennot model of lozenge tilings. (Not all 
orientations are shown.) 



Let Y(X/ll\x) be the graph exemplified in Figure 17 It 
is a diamond grid of height a and width X { + b. All edges 
point down, and those in the ith row that point southwest have 
weight x,. On the top row there is a left endpoint at each po- 
sition b + i — \i\ and on the bottom a right endpoint at each 
position b + i — X[. Note that the left endpoints are at all of the 
positions not of the form \i x ■ + b + 1 — i and likewise the right 
endpoints are at the positions not of the form One can check 
that the Gessel-Viennot matrix of X(X/jJ.;x) is D(X/fi;x). At 
the same time, noting the gaps between the marked endpoints, 
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Figure 18: An Aztec diamond of order 3 



Domino tilings of an Aztec diamond are a well-known ana- 
logue of lozenge tilings of a hexagon. Recall that an Aztec 
diamond of order n is the polyomino consisting of those unit 
squares lying entirely inside the region \x\ + \y\ < n + 1, as 
shown in Figure A domino tiling of an Aztec diamond 
corresponds to a matching of the graph dual to the tiling by 
squares; an example of such a graph is the one on the left 
in Figure [|. Denote this graph by Z A (n) and let M A {n) be a 
Kasteleyn-Percus matrix for it. 

Theorem 16. If M A (n) is a Kasteleyn-Percus matrix for 
domino tilings of an Aztec diamond of order n, then 

cokerM A (n) = Z/2 © Z/4 © ■ ■ ■ © Z/2". 

Theorem [l6| extends the result that the number of domino 
tilings of an Aztec diamond of order n is 2"'" +1 )/ 2 \ 



Proof. We give two arguments, one using Kasteleyn-Percus 
matrices and the other using Gessel-Viennot matrices. First, 
define a binomial coefficient matrix B(n) by 



B (")o<U< 



12 



(We will assume that rows and columns of other matrices are 
numbered from as well.) Define annxn+1 matrix L(n) by 
putting the n x n identity matrix to the left of a null column, 
and define R{n) by putting the nxn identity matrix to the right 
of a null column. For example here are B(4), L(3), and R(3): 



(l 0\ 
110 
12 10 
\l 3 3 1/ 

One can check that the matrix 



M A {n) = R{n) ® R(n) J + L{n) <g> R(n) T 

+ R{n)®L(n) T ' -L(n) ®L(n) T 

is a Kasteleyn-Percus matrix for the graph Z A (n). At the same 
time, 
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B(n)L(n)B(n+l) 



L{n) 



and 



B(n)R(n)B(n + iy l =R{n)-L(n). 
It follows that 

M'(n) = (B(n) ®B{n + \)- T )M(n)(B{n + l)" 1 ®B(n) T ) 
= R(n) ®R(n) T - 2L(n) ®L{nf . 

Since B(n) is a triangular matrix, M'(n) is equivalent to M(n). 
On the other hand, one can check that, modulo permuting 
rows and columns, 

M\n)=X{\)®X{2)®---®X{n) 

©F(i)©y(2)®---©y(n), 



where 
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Figure 19: Gessel-Viennot model of an Aztec diamond. (Not all 
orientations are shown.) 



the theorem follows. 

In the Gessel-Viennot approach we first observe that if G(n) 
is a graph of the type in Figure [l9]then its transit-free resolu- 
tion G'(n) as described in Theorem ^ is the Aztec diamond 
graph Z A (n). (The left and right endpoints of G(n) are num- 
bered in the figure, with the left endpoints on top. The first 
left and right endpoints coincide, which is degenerate but al- 
lowed.) Let V(n) be the Gessel-Viennot matrix of G(n). The 
entry V(n) { ,■ is the Delannoy number D{i,j) j2^ , §6.3] (see 
also Sachs and Zernitz [p4{]), because G(n) matches the defin- 
ing recurrence 

D(i, j) = D(i, j - 1) +D(i - 1 , j) +D(t - 1 , j - 1) 

D(0,i)=D(i,0) = 1. 

A standard formula for Delannoy numbers is 

o(u)=E' 



In matrix form this identity can be expressed 

V(n) =B(n)V'(n)B{n) T , 
where V'(n) is the Smith normal form of V(n) with 



Thus the Smith normal form of V{ri) also establishes the the- 
orem. □ 

Finally, we could at least conjecturally extend Theorem p"o| 
with the same variations as those we considered for lozenge 
tilings: Domino tilings with symmetry, ^-enumerations, im- 
possible enumerations, Aztec diamonds with teeth missing, 
etc. For example, Tokuyama established a relation be- 
tween generating functions of lozenge-type and Aztec-type 
Gelfand triangles. We conjecture that this identity can be 
extended to an equivalence of Kasteleyn-Percus (or Gessel- 
Viennot) matrices. We leave this and other possibilities to fu- 
ture work. 



Since 



cokerX(k) = cokerT(fe) = Z/2* 5 



6. APPENDIX: SMITH NORMAL FORM 
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Theorem 17. If M is a k x n matrix over a principal ideal 
domain R, then there exist invertible matrices A and B such 
that 



If we pre-multiply by 



Sm(M) =AMB 

is diagonal and Sm(M) ; ; divides Sm(M) (+1 . , j. 
Sm(M) is a Smith normal form ofM. 

Evidently 
cokerM = cokerSm(M) 



The matrix 



■■R k - n ®R/Sm(M) 11 < 



)R/Sm(M) n>n 



if k > n and 



cokerM = cokerSm(M) 
= R/$m(M) 11 Z 



)fl/Sm(M) 



k.k 



otherwise. (If R = Z and we allow n, but not £, to be infinite, 
then Theorem |l7| is equivalent to the classification of finitely 
generated abelian groups.) It is easy to show that this decom- 
position of cokerM is unique and that Sm(M) is unique up to 
multiplying the entries by units. 

Proof. We argue by induction on n or k. Since R is a PID, any 
two elements a and b have a greatest common divisor 

gcd(a,b) — ax + by 

which is unique up to a unit factor. Moreover, R is Noethe- 
rian, which means that any chain of divisibilities a j+l \a t must 
eventually be constant up to unit factors. Thus we can argue 
by induction with respect to divisibility of non-zero elements 
of R. Since we can replace M by any equivalent form, we as- 
sume that the entry M l l ~ a is not divisible by any entry in 
any form of M, except for those that are a times a unit. We 
claim that a divides every entry of M. Otherwise there is an 
entry M ; • = b such that a j b and b \ a. If £> is in the first row or 
column, then after permuting rows and columns and possibly 
transposing M, M 



1.2 



M : 



Let c = ax + by be a greatest common divisor of a and b. If 
we post-multiply M by the matrix 





fx 


a/c 





B = 


y 


-b/c 




K 








the result is a form of M with the entry c, which contradicts 
the choice of a. If a divides every entry in the first row and 
column but not some entry b elsewhere, then we first pivot at 
the position (1,1) to obtain: 



M 



( a 
b 



\ 





f\ 1 




A = 


1 






V o 





and post-multiply by B above, we again produce the entry c. 

Given that M l l — a does divide the rest of M we perform a 
deleted pivot at ( 1 , 1 ) and assume normal form for the remain- 
ing submatrix by induction. □ 



Theorem 18. If A is an alternating nxn matrix over a prin- 
cipal ideal domain R, then there exists an invertible matrix B 
such that 

Sm„(A) = B T AB 

is block-diagonal with 2x2 blocks, and such that 
Sm„(A) 2; 2i+\ divides Sm a (A) 2i - +2 2H-3- The matrix Sm a (A) is 
the alternating Smith normal form of A. 

Proof. The argument is the same as the one for Theorem [l7| 
We assume that A l 2 = a is minimal among all entries of all 
forms of A. If it does not divide some entry in the first two 
rows and columns, then after permuting rows and columns, 
that entry is A l 3 = b: 

( a b 
-a • 
l= ~b ■ 



Let c = ax + by be a common divisor and let 



B = 



( 1 

Ox a/c 
y -b/c 



V 



\ 





7/ 



Then B T AB has the entry c, a contradiction. If a divides every 
entry in the first two rows and columns but not some other 
entry b, then we c an perform a symmetric pivot at ( 1 , 2 ) . After 
permuting rows and columns A then has the form: 



A': 



/ a 
-a 


-fc 



V 



Let c — ax + by be a common divisor and let 



B 



( 1 
10 
10 10 
1 

V o 



o 
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Then B AB has the form of the previous case. form for the remaining submatrix. □ 

Finally if a divides every entry of A, we perform a deleted 
symmetric pivot at (1,2) and inductively assume the normal 
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